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Abstract. We use a connection between relativistic hydrodynamics and scalar field theo- 
O ! ry to generate exact analytic solutions describing non-stationary inhomogeneous flows of 

2P ' the perfect fluid with one-parametric equation of state (EOS) p = p{e). For linear EOS 

p = KE we obtain self-similar solutions in the case of plane, cylindrical and spherical sym- 
' metries. In the case of extremely stiff EOS (k = 1) we obtain "monopole + dipole" and 

"monopole -I- quadrupole" axially symmetric solutions. We also found some nonlinear EOSs 
that admit analytic solutions. 
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^ ; 1 Introduction 

> , 

, Relativistic hydrodynamics (RHD) is extensively used to describe various processes from mic- 

I roscopic to cosmological scales. Astrophysical applications of RHD deal with powerful non- 

stationary phenomena such as the gamma-ray bursts and hypernovae explosions. On the other 
Q>^ . hand, RHD solutions are sources of hydrodynamical models of elementary particle multiple 

I production. Though at present there are powerful methods to find numerical solutions of hydro- 

dynamical problems, the price for exact analytic representations of RHD flows is still high. 
However, the nonlinear structure of the RHD equations hampers the quest of analytic solutions 
and this is the reason why RHD problems having exact solutions typically deal with simple 
^ ' equation of state (EOS) such as p = ke (p is the pressure, e is the proper frame energy density, 

K = const). This linear equation of statqj is often used in multiple pion production theory since 
works of Landau [1] and Khalatnikov [2] . 

The most part of known exact RHD solutions concerns plane (1 + 1 dimensional) barotropic 
relativistic flows. The Riemann simple waves are the example (see, e.g., [3j). In the case of the 
linear EOS the method of Khalatnikov [2] allows to obtain a general solution in an implicit form. 
However, the solutions obtained by means of this method have rather complicated representation. 
This motivated different authors to look for more simple solutions to be used in models of 
multiple pion production O El [71 ; these solutions have extensions to the case of spherical and 
cylindric symmetry. Recent interest to such solutions is due to high energy heavy ion collisions 
(see, e.g., PHOllIl]). 

The three-dimensional case of RHD is less developed. No general analytic solution is known 
in the case of relativistic spherically symmetric flows as well as in the case of cylindric sym- 
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'^In the case of the EOS p = ne+B the constant B can be easily eliminated from RHD equations by redefining e. 
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metry. The best known in theory of astrophysical relativistic outbursts is the Blandford & 
McKee approximate self-similar solution [T2] that describes ultra-relativistic blast waves of the 
fluid with the linear EOS; for modifications and generalizations (also in the ultra-relativistic 
approximation) see [131 E] • 

Various techniques were invented in classical hydrodynamics in order to find analytic solu- 
tions. In particular, one may consider potential flows of ideal fluid to reduce the number of 
unknown functions. In this paper we study RHD flows that may be considered as a relativistic 
analog of the classical potential flow: in these flows the four-velocity is proportional to a gradient 
of a scalar function. The relation between RHD flows and the scalar fields is well known (see, 
e.g., [HEllll]) and at present it may be considered as a part of physical folklore. This relation is 
extensively used in cosmology in connection with the dark energy problem |16l [T71 \TE[ [T9] . Here 
a number of unusual equations of state has appeared; they have been considered in the context 
of spatially homogeneous cosmological systems (see, e.g., [TO] for the review). 

In this paper we use the connection between RHD and the scalar field theory that allows 
us to reduce the problem to an equation for a single scalar function (instead of four-velocity 
and energy density). The aim of this paper is to study the ability of this trick to generate 
analytic RHD solutions that represent non-stationary and non-homogeneous flows of relativistic 
perfect fluid. The paper is organized as follows. In Section [2] we present basic equations 
and relations between the scalar field and hydrodynamical variables. In Section [3] we present 
solutions for an extremely stiff EOS: there are plane solutions, spherical solutions including 
different combination of outgoing and ingoing spherical waves; at last we present examples of 
three-dimensional solutions that are obtained from spherical harmonics of the scalar field. In 
Section H] we present some nonlinear equations of state that allow us to find exact RHD solutions. 
In Section [5] we present two families of self-similar solutions in the case of the linear EOS. In 
summary we briefly discuss the results obtained, their relation to the other known solutions and 
possible generalizations. 



2 Notations and basic equations 

The relativistic equations of the perfect fluid dynamics may be written as conservation laws 

(e.g., my- 

duT'^" = 0, (1) 
where 

T^- = (p + e)n'^n'^_p(^M-, n'^n^ = l, = diag (1, -1, -1, -1) ; (2) 

u^^ is the four- velocity, p is the pressure and e is the proper frame energy density, the Greek 
indices run from to 3, the speed of light c = 1. In the general case equation ([T]) must be 
supplemented with the baryon number conservation equations. However, it is not necessary in 
this paper because we deal only with the one-parametric EOS p = p(e); provided EOS be known, 
the system ([TJ is complete. Further we concretize different equations of state in every section. 

If we have a conserved energy-momentum tensor of any field with known solutions, we can 
use it to find a hydrodynamic solutions for e, u'^. However, to do this one needs to represent 
this energy-momentum tensor in the form ([2]); moreover, in order that these solutions represent 
physically admissible motion of the relativistic fluid one must provide the time-like behavior of 
the four-velocity. This is not always possible. However, solutions with a spacelike are also 
of some interest: they demonstrate that formal solutions of completely relativistic dynamical 
equations may admit tachyonic motions. 
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Transition to the proper frame of a local fluid element shows that arbitrary tensor can be 
represented in the form if the matrix has only one eigenvalue e > corresponding to 
a timelike eigenvector and triple-degenerate eigenvalue p with space-like eigenvectors. It is not 
always possible to fulfill these requirements for an arbitrary Tlf. However in the case of a scalar 
field with the Lagrangian 

L = F{S), 5=^(^,,(^," (3) 
the corresponding energy-momentum tensor 

Tap = ^,aV,pF'{S) - ga/3F{S) 

takes on the form (l2|), if we put [J: 

p = F{S), e = 2SF\S)-F{S), (4) 
Ua = ±ip,a{2Sr^/\ (5) 

The equations ([4]) provide an effective EOS that is represented parametrically. 
The equation of the scalar field following from ([3]) is 

=0. (6) 

Formulae ([H), ([5|) represent an admissible hydro dynamic flow provided that 

e>0, S>0. (7) 

This is the only condition to check, if we have a solution of equation ([6]). 

3 Extremely stiff EOS: plane and spherical solutions 
3.1 Plane solutions 

The most simple is the case of a massless scalar field corresponding to the extremely stiff EOS: 

F{S) = S, p = e = S. 

This EOS is considered throughout the whole section. In this case we deal with a linear wave 
equation for the scalar field: 

df.d'^ip = 0. (8) 

In the case of plane symmetry [x^ = t, = x) it is easy to see that the general solution 
ip = ip {t — x) + xit + admits a hydrodynamics interpretation if only 

^l^' {t - x) x' {t + x) > 0; (9) 

this makes it evident that the solutions corresponding to the only simple wave moving in one 
direction have no hydro dynamical counterpart. For hydrodynamical variables we have 

dipldx lb' (t — x) — V it + x) ,,, . s . X 

^ = — = -7rW = T777 L ;L ' e = 2i,'{t-x)x'{t + x). 10 

u'^ oip/ot ip' [t — X) + X [t + X) 

For -0' (x) = x' (x) = Ax~^ (throughout the paper A = const G M, ^ 7^ 0) we obtain the 
scaling solutions [3 El [3 [8] : v = x/t, e = 2A'^(t'^ - x'^)~^ that are well defined inside the light 
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cone: < t. For x > t > the hydro dynamical interpretation fails. If the solution for t = ti, 
|x| < xi < ti is considered as initial data for further hydrodynamic evolution, either the solution 
must be complemented by correct hydrodynamic data for \x\ > xi, or the equations must be 
modified (cf., e.g., [HI [22]) to extend the solution for all x. Note also that after rescaling A ^ iA 
we obtain a formal "tachyonic" solution outside the light cone with e > 0, but with v > 1. 

Choice ^' (x) = x' (^) = Ax^, \x\ < t (here and below A G M, A 7^ 0), generates similar 
solutions with the energy density e = 2A'^{t^ — x^)"^ and \v\ < 1. For A < we have an outflow 
{v > 0); this solution is singular at the light cone. For positive A we have an inflow (u < 0); in 
particular, for positive even A the energy density has no light cone singularity. Simple rescaling 
of factors in the functions tp, x generates physical or tachyonic flows correspondingly inside or 
outside the light cone. 

In order to construct a physical flow outside the light cone one can choose x' i^) = ^' = 
Ax^, X > t > 0. We have an outflow for A < 0, inflow for A > 0. This case can be easily 
extended to negative x and/or to negative t. For A = 1 we obtain e = 2^^(x^ — t^), v = —t/x. 
This is a kind of external scaling solutions discussed in the papers [Oj [10]. In order to obtain 
strict inequality Q and regular solutions Va;, t, the above power-law choice for tp and x niay 
be replaced, e.g., by x' i^) = 4'' ~ exp(7x"'), where 7 = const S M, n is a positive integer. 

3.2 Spherical solutions 

The general solution in the case of spherically symmetry (x" = t, x^ = r) is 

V = r'^ [ij{t-r) + x{t + r)], 

where and x arbitrary functions. Here, as distinct from the planar case, it is possible to 
use the outgoing wave solutions (with x = 0): 

V= , (11) 



1,2 r 



P = ''~ 2r4 



2rih' 
1+ ^ 



whence 

I H 

In order to provide ([7|), it is necessary 

lb' (a) , , 

2^^r<-l, a = t-r, (12) 

therefore the function |^ (a) | is decreasing. Evidently it is impossible the condition ()12p to 
hold for all t, r (in particular, it is violated for r ^ 0) and the hydrodynamical interpretation 
of T^jy corresponding to the solution (jlip is possible only in a bounded region. When the sign 
of S changes and e = this solution must be matched to vacuum in the way appropriate to 
hydrodynamic flow (see Appendix A). 

Example. Consider a solution 

Ci > and C2 > being real constants, n is positive integer. In this case ([7|) is fulfilled if 

This holds at least in the domain {(t,r) : t > r > t/ (2n + 1)}. As we see from Appendix A 
this solution is matched to vacuum for t — r = ti, if Citj"" = C2. For r > t the hydrodynamical 
interpretation fails and it is either necessary to match the solution through a discontinuity, or 
the solution is destroyed by an external perturbation having trajectory r = —t + const. 
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3.3 Combination of outgoing and ingoing spherical waves 

In this case, in order to provide regularity of the solution for r — > 0, we put x (0 = (*) 

if [t, r) = {i/j {t - r) - tp {t + r)} . 

Consider the power-law choice of ip (x) = —Ax~'^, A > 0, n > 0. For every n the inequalities d?]) 
must be analyzed separately. Below we present the hydrodynamic solutions corresponding to 
four natural values of n; they have a hydrodynamic interpretation inside the light cone t > r. 
For n = 1, t > r: 

¥^ = -^2372' '=P=J^^^ ^ = -f (13) 

Appropriate rescaling of the constant A yields a tachyonic solution outside the light cone. 
For n = 2, t > r: 

4tA _ _ 8^2(9^2 _ ^2) ^ 

'^-~(t2_^2)2' ^-P- (^2_^2)5 ' ^ - 3^2 + ^2 • (^^) 

For n = 3: t > r: 

+r2 _ _ 32^2(9*^ + 2*^2 + r^) _ r(5*2 + r^) 

^--2^(t2_^2)3' ^-P- (t2 _ r2)7 ' ^ - 3t(t2 + ^2) " (l^) 

For n = 4, * > r: 

= -8^(i2_,2)4' 



_ 32^2[;^o^4(^2 ^ ^2) ^ ;^5^2(^4 ^ ^2^2 ^ ^4) _ ^6] 2tr(3r2 + 5*2) 

^-•P- (*2 - r2)9 ' 5t4 + io*2r2 + r4- ^ ^ 

In a more general case {n > 1, yl > 0) one can show that hydrodynamic flow is correctly 
defined at least in the domain {(*,r) : k„* < r < *}, where = {an — l)/(on + l)j o,n = 
[2{n + 1)]^''"', and in this region > 0. 



4 Extremely stiff EOS: solutions based on spherical harmonics 

In this section we proceed with the extremely stiff EOS; here generation of RHD solutions 
utilizes the well-known formula for the solution of the wave equation ([8]) using spherical functions 
Yim (^) f) (in spherical coordinates) 

oo I 

^{t,r,9,ip) = Y, E fim{t,r)Yi„,{e,ip), (17) 

1=0 m=~l 

Jlm = r -— < 

\r or / \ r 

where (* — r), (* + r) are arbitrary functions. 

Investigation of inequality d?]) for (jl7p is rather complicated in the general case. However, 
new special solutions obeying d?]) may be easily obtained as follows. Take one of spherically 
symmetrical solutions p^ - p^ from the previous section and add some items from the sum 
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in the r.h.s. of (jl7p . These items must be sufficiently small in comparison with (jl3p - (|16p to 
preserve inequality ([T]) at least within a compact domain inside the light cone. This is possible 
for a special choice of appropriate functions (jlSp . Below we present two such special solutions 
that represent a hydrodynamical flow within the light cone t > r. 

In the case of axial symmetry we obtain the three-dimensional velocity components from the 
relations 



The next case corresponds to the monopole + dipole contribution into ip: 

,t — br cos 9 
ip = A- 



(f - r2)2 



6 is a real constant. If |6| < 1, then the condition d?]) is satisfied for all t > r. Corresponding 
solution in terms of hydrodynamic variables is 



A^ 



P 



2{f _ ^2^5 



{l-h^){f -r'^) + S {t-hr COS of > 0, (19) 



^ _ 4tr- (3r2 +t2)fecos(g) q _ {t^ - r'^)b sm{9) 

^ ~ 3*2 + r2 -4tr6cos(6l) ' " ~ r[3t2 + r2 - Atrb cos{e)]' ^ ^ 

The other solution is generated by a monopole + quadrupole contribution 

3*2 + r2 + 6r2p2 



If = A- 



(t2 



where P2 = [3cos2(0) — l]/2, 6 is a real constant. We checked that the condition ([7]) is also 
satisfied in this case for all t > r, at least if |5| < 1. Corresponding hydrodynamical solution is 

e = P= 2^^2^^^2y [I44i' + 16r^(2t2 + r^) (i + 26P2) 

- 2rH\P2 + - r2) + m'^P^r^] > 0, (21) 
r 10*2 + 2r2 + (t2 + 2r2)6P2 e b {t^ - r^) sm{e) cos{e) 



St 2*2 + 2r2 + 6r2P2 ' 2t 2t^ + Ir"^ + br'^ P2 

Nonlinear barotropic EOS 



(22) 



In this section we consider equation ([6j) in the case of plane (n = 0), cylindrical (n = 1) and 
spherical (n = 2) symmetry. The dependence p = p{e) will be specified later. In this case 
equation ([6]) can be written as 

We are looking for solutions of the equation ([6]) of the form 

^ = ^{a), a={e-r^)/2 (24) 
In this case 
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and substitution into equation (|23|) yields 



d 

2a — 

da 



da 



+ {2 + n)F'{S)^ = 0. 

da 



(25) 



It is convenient to introduce a new variable r by the relation a = S = {l/2){dip/dT)'^ . 

Then equation (p5]) yields 



t1+"F'(5) 



dip 



const. 



(26) 



We investigate the most simple cases, when equation ()26p can be easily solved with respect 
to dip/dr and the EOS is expressible in terms of elementary functions. 
Consider first the EOS 



p = e 



ln{ — ]+ B 



F{x^/2) = eQx[\u{x) + B], 



(27) 



where i? is a dimensionless constant, £q is a constant having the dimension of the energy density. 
Substitution into equation ([26]) gives us 



dip 



C 

-l+n 



B-l], 



C being an arbitrary real constant. The solution of hydrodynamical equations with EOS (|27|) is 



v = r/t, e = eoexp I - S - 1 



C 



Vt2 



t> r. 



(28) 



This solution is formally valid for interior of the future light cone. However, for usual hydrody- 
namical interpretation we neecH that 



cl = dp/de £ [0, 1], 



(29) 



where Cs is the speed of sound. This yields an additional restriction of the domain of RHD 
solution: e'^ < ee^/eo < 1 whence t"^ - r'^ > C2/(i+"), C > 0. 
In a more complicated case we consider 



N 



F{S) = eoVS 



m=0 



{-1)"'N\ 

(iV-m)! 



lnV5 



N-m 



with 



dF 



N 



dS 2^/S 

Corresponding EOS is represented parametrically {R = ^/S) as 



N 



N 



(-l)"^iV! 
p = £qR 2^ — — — {InR} 



and both e and p may be expressed as functions of enthalpy 
e + p = £QR{\nR)^ . 



(30) 



^See, however, remarks in the last section. 
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In this case 

2 dp hi R 



(31) 



' de N 
On account of equation (j26p 

t1+" {\nRf = const. 
This yields the hydrodynamical flow for the EOS (j30p with 

v = r/t, i? = exp [c(t2-r2)-(i+")/(2^)j , < r < t, (32) 
C is a real constant. Taking into account of equations (j31|) . (|29|) yields 

1 < < e^, t2 _ ^2 > (^/^)2JV/(l+n)^ ^ > Q_ 

6 Linear EOS p = ne 

It is easy to find the Lagrangian ([3]) corresponding to the EOS p = Ate, where k = c^, < k < 1. 
Using ([3D we obtain a simple differential equation for F yielding 

F{S) = 5", a = (33) 

2k; 

up to inessential multiplier. 

Such Lagrangians and their generalizations have been discussed in [6l [T6t [T71 123| . Now we 
consider equation Q in the case of plane, cylindrical and spherical symmetry on account of 
equation ([33l) . Equation (p3|) is reduced to the following equation 



First, we are looking for the solutions of the form (j24p . We have a special case of equation (j26p 
that can be easily solved yielding the famous scaling solutions of relativistic hydrodynamics 

13 El [3 i: 

c 

^=—fTTUuT^^ v = r/t, 0<r<t, C > 0. 

^(l+fe){l+n) ' ' ' ' 

Now we shall look for solutions of equation (|34p of the form if = i^i^,), £, = r/t. Then 

^-^^[d^j ■ 

One can check by direct calculations that equation (j34p leads to the equation 



dz nz 



di ae-iy 

For ^ > 1 we obtain 
dip C 



C = const. 



2a-l 
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r 



Figure 1. Energy density e{t, r) according to equation ((35|) . rt = 2, for different times t = ti in the case 
of fcn < 1 (lines 1, 2: ti > and kn > 1 (lines 3, 4: ^3 < t^). 

The hydrodynamic variables are as follows 
t 

^ = p = Ke = S°', S = 2^2nK(^2 _ i2y-nK'^ ^^^^ 

the solution is valid outside the light cone (0 < t < r). Again, after rescaling of the constant C 
we have a tachyonic solution inside the light cone with e > 0, \v\ > 1. 

For n = 1,2 there is the light-cone singularity for all k G (0,1). The same situation occurs 
in the case of k < 1/2 for n = 2. In the case of a continuous flow the solution (|35|) may be 
matched with a numerical solution through a sound wave. 

In the case of spherical symmetry (n = 2) and fc > 1/2 the energy density has no light cone 
singularity and we have a continuous solution (j35p that is matched with vacuum at the light 
cone (Fig. [T]); this may be checked using equation (j36p from the Appendix El In the case of 
kn = 1 the density does not depend on time describing an outflow supported by an energy 
source from the interior. 

7 Summary and discussion 

In this paper we considered relativistic potential flows of the perfect fluid that permit reductions 
to an equation for one scalar field. The restriction to the potential flow allowed us to derive 
analytic RHD solutions that are interesting from theoretical viewpoint, e.g., in the theory of 
elementary particle multiple production; they may be used to test numerical RHD computer 
routines. The relation between RHD and the scalar field proved to be most fruitful in the case 
of the extremely stiff EOS. In this case the problem is reduced to the linear wave equation 
for the scalar field yielding hydrodynamical solutions in the cases of planar, cylindrical and 
spherical symmetry. Moreover, a family of "non-symmetric" RHD solutions may be generated 
using multipole expansion for the scalar field; this is illustrated by examples that arise from 
combination of lower multipoles. We also found a family of nonlinear EOSs that permit analytic 
solutions defined inside the light cone and having behavior like that of famous scaling solutions 
[6l O [8] . In the case of the linear EOS p = ke some of our solutions reduce to these solutions 
O El [71 E] . Also, we derived the solutions outside the light cone from self-similar solutions of 
the nonlinear scalar field equation. 

Most of our solutions have light cone singularities like that of O El O [8] . These solutions 
cannot be directly matched with vacuum and may be used to describe a flow in a finite volume 
until it is destroyed by perturbations moving from the boundaries. Nevertheless, there are 
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exceptions when regular solutions in a whole space may be generated, especially in the case of 
k = 1 and in the case of spherical solutions (j35|) with k > 1/2. 

After we have sent the first version of our manuscript to the journal the paper [lOj was 
publishecH containing new interesting RHD solutions. Some of the solutions obtained in the 
present paper were independently found in [10] using different methods. In particular, our 
solution (|35p corresponds to equations (24), (25) of [10] that have been obtained as an extension 
of "accelerating" solutions of [9]. Also, our solution (jlOp in the case of the plane flow and 
extremely stiff EOS is the same as solution (114), (115) of |10j . 

To outline possible generalizations, we note that the condition about the speed of sound 
(0 < A; < 1 in the case of the linear EOS or equation (|29p of Section [5|) may be relaxed. This 
may be of some physical interest: the EOS with negative pressure, negative k < or k > 1 are 
discussed in cosmology [13 [19], see also [21] concerning the possibility of "superluminal" EOS. 
The results of the Section [5] will be extended, if we do not restrict ourselves to simplest EOS in 
terms of elementary functions. In Section [6] we confine ourselves to the self-similar solutions of 
the field equation ()34p of the form ip = ip(r/t) thus yielding an ordinary differential equation. 
However, more general substitution ip = r°'ip{r/t) also leads to an ordinary differential equation, 
though more complicated. 



A Matching of solutions with vacuum 

To perform this matching in a hydrodynamical way we use the condition of zero energy- 
momentum flux through a boundary with vacuum 



T^k, = 0, (36) 
where the vector ky is orthogonal to the boundary hypersurface. This yields 

Simple analysis shows that this is possible only if at the boundary 



^,a'P'^ = 0, '■P,ak'^ = 0. 

In the two-dimensional case we have ip\ = ip\. 

Consider the light-like surface t — r — ti = 0, q = ti, then k^ = k^ = —ki = 1. In the case of 
the solution (|lip the boundary condition is 



that is at the boundary with vacuum ip {ti) = 0. 
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